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ABSTRACT 

The recently detected 7-ray emission from Starburst galaxies is most commonly considered 
to be diffuse emission arising from strong interactions of accelerated cosmic rays. [Mannheim| 



et al. (2012 1, however, have argued that a population of individual pulsar-wind nebulae 
(PWNe) could be responsible for the detected TeV emission. Here we show that the Star- 
burst environment plays a critical role in the TeV emission from Starburst PWN, and perform 
the first detailed calculations for this scenario. Our approach is based on the measured star- 
formation rates in the Starburst nuclei of NGC 253 and M 82, assumed pulsar birth periods and 
a simple model for the injection of non-thermal particles. The two-zone model applied here 
takes into account the high far-infrared radiation field, and different densities and magnetic 
fields in the PWN and the Starburst regions, as well as particle escape. We confirm that PWN 
can make a significant contribution to the TeV fluxes, provided that the injection spectrum of 
particles is rather hard and that the average pulsar birth period is rather short (~ 35 ms). The 
PWN contribution should lead to a distinct spectral feature which can be probed by future 
instruments such as CTA. 
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1 INTRODUCTION 

Starburst (SB) galaxies have a very high star-formation rate (SFR) 
and hence supernova (SN) rate in a very localised region - the SB 
region - which is often located in the centre of the galaxy. The 
paradigm of cosmic-ray acceleration in supernova remnant shock 
fronts and the very large amounts of gas in the SB region, avail- 
able for proton-proton interactions and subsequent 7r° -decay 7- 
ray emission, made them suspected emitters of high-energy (HE; 
100 MeV < E sC 100 GeV) and very-high-energy (VHE; E ^ 
100 GeV) 7-rays(e.g. fV61k et al.|1989)|Akyuz et al.|1991) emis- 
sion. Indeed, 7-ray emission from the prototypical SB galaxies 
NGC 253 and M 82 using the space-based Fermi-LAT ( |Abdo et al.] 
|2010[ > and the ground-based Cherenkov telescopes H.E.S.S. and 
VERITAS ( |Acero et al.|2009l| VERITAS Collaboration|2009| > has 
recently been reported. There are no indications of 7-ray variabil- 
ity, and in the one case where the instrumental resolution is suffi- 
cient to discriminate, the TeV signal from NGC 253, the emission 
is compatible with originating in the SB nucleus, without a signif- 
icant contribution from the disc of the galaxy ( Abramowski et al. 
|2012l 

Interpretations of the detected GeV and TeV 7-ray emission 
from SB galaxies are possible within the traditional framework if 
~ 10% of the kinetic energy of SN explosions is converted into 
particle acceleration and about (20 — 30)% of these protons lose 
their energy in interactions with target material (see e.g. |Lacki| 



l et al.||2011[|Ohm & Hinton||2012| |Paglione & Abrahams||2012{ 
Abramows ki et al.||2012| >. Whilst this scenario is certainly plau- 
sible, it may be somewhat oversimplified. For example, it is not 
clear that a contribution from individual sources can safely be ne- 
glected at all energies. The TeV 7-ray source population observed 
in our own Galaxy is diverse and shows emission from 7-ray bina- 
ries, stellar clusters, pulsar wind nebulae (PWNe), as well as SNR, 
with PWNe being the dominant object c lass (e.g. |Hinton & Hof-| 
|mann| [2009 ). Based on this observation, Mann heim et al.| ( 2012) 
argued that a major part of the emission seen in the TeV band from 
NGC 253 and M 82 could be explained by the combined emission 
from individual PWNe. The authors used the Galactic population of 
TeV-detected PWNe to estimate the contribution of leptonic 7-ray 
emission to the signal seen from both SB galaxies. However, the 
environment in which most of these PWNe evolve is significantly 
different from the physical conditions in SB nuclei. In this work 
we follow the idea of Mannheim et al.|{ |2012| >, but study the PWN 
population and the resulting non-thermal emission in the SB based 
on the SFR and time-dependent pulsar injection power, rather than 
the (incomplete) sample of Galactic TeV PWNe. 

2 PWN IN A STARBURST ENVIRONMENT 

Pulsar wind nebulae are formed when the cold ultra-relativistic 
wind from a pulsar interacts with its environment. Particles in 
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PWNe are accelerated up to ~ 1 PeV and produce synchrotron and 
Inverse Compton (IC) emission from radio to VHE 7-ray energies 
in interactions with magnetic and radiation fields, respectively (e.g. 
|Gaensler & Slane|2006| >. The majority of PWNe in the Milky Way 
evolve in the Galactic disc, where the target radiation fields for the 
IC process are the Cosmic Microwave Background (CMB), visible 
starlight, and reprocessed starlight in the infrared. The typical en- 
ergy densities of these radiation fields are U p h ~ leVcm" 3 . The 
dominant radiation field in the nuclear region of SB galaxies, on the 
other hand, is far-infrared (FIR) radiation with orders of magnitude 
higher energy density. Melo et al. (2002) for example found a total 
FIR luminosity of 2 x 10 iu L Q in the SB region of NGC 253. The 
NGC 253 SB has a cylindrical shape with a radius of ~ 150 pc and 
a half-height of ~ 60 pc i jWeaver et al.|2002| ). Assuming that each 
point in the cylinder radiates with the same emissivity at the tem- 
perature of the cold dust (^50 K), £/ p h,NGC253 — 2400 eV cm" 3 , 
a factor 2.9 higher than for a point source with the same luminosity 
at 150pc distance]^] This value is more than three orders of mag- 
nitude higher than is typical in the ISM and implies that the cool- 
ing time of highly-relativistic electrons in the SB region is much 
shorter, i.e. ~ 200 yrs compared to ~ 3 x 10 5 yrs for a 1 TeV elec- 
tron. IC cooling in such an environment will dominate over syn- 
chrotron losses unless B > 500 fiG. 

The central source that powers the PWN is a pulsar that spins 
down and converts a fraction of its rotational energy into non- 
thermal particles and subsequently into radiation. Over time the 
pulsar's rotation slows down, and the energy put into the nebula 
decreases. The spin-down luminosity evolves in time as 



E(t) = E /(l + t/r) p , with r = P /P (n - 1), 



(1) 



where r is the characteristic spin-down time, given by the birth 
period of the pulsar Po and the first derivative of Po. p = (n + 
l)/(n — 1), where n is the braking index and is measured only for 
a handful of pulsars, where it lies between ~ 2 and 3. In order to 
investigate the effect of the different radiation fields we will apply a 
single-zone, time-dependent model, where particles are injected ac- 
c ording to Equation[T] T his approach is similar to the one followed 
in |Hinton & Aharonian|(2007) . 

Figure[TJshows the fraction of the total injected electron power 
that is radiated in different spectral bands for a PWN evolving in 
both SB-like, and ISM-like environments for different ages. Parti- 
cles are injected according to a power-law in energy N(E) oc E~ a 
with index a — 2.0. Recent studies of prominent Galactic PWNe 
(Zhang et al. 2008] e.g.) suggest that the conversion efficiency of 
spin-down power into non-thermal particles e is ~ 1/3, with a re- 
cent lower limit of 0.3 derived for MSH 1 5-52 |Schock"etlLl ( [2(JT0l ). 
However, the fraction of spin-down power that is converted into 
magnetic fields in young PWN is known to be small (Kennel & 
Coroniti| 19841 |Del Zanna et al.|2006||Gelfand et al.|20091|Tanaka| 



& Takahara|2010 



e.g.), so that e may approach 1. Non-radiative 
losses and/or complex spectral shapes for the lower energy parti- 
cles may be responsible for the apparently lower injected power in 
particles (Spitkovsky 2008 , Slane et al. 2010, e.g.). In the following 
we assume a constant value of e = 0.5 and neglect non-radiative 
losses, with an estimated uncertainty of a factor of 2 in either direc- 
tion. 

Figure[T]illustrates that emission from an object evolving in a 
typical ISM region is likely to be synchrotron-dominated at early 



times, when the B-field is high. At later stages, when the mag- 
netic field decreases as the nebula expands, the fraction of energy 
converted into IC emission is expected to increase. This effect pro- 
vides a plausible explanation for the number of bright, middle-aged 
PWNe detected in VHE 7 rays that have rather low X-ray flux (see 
e.g. |de Jager & Djannati-Atai'|2009| >. It has been suggested that such 
(more numerous) middle-age PWN may be the main contributors 
to the VHE emission of SB galaxies (Mann heim et al.|2012[ >. How- 
ever, the energy output of a Crab-like pulsar after 10 3 years is ex- 
pected to be a factor of ~ 50 higher than at an age of 10 4 years, 
more than compensating for the smaller number of young objects. 
In a typical ISM environment these younger objects would be ex- 
pected to be synchrotron-dominated, with low IC efficiency. Due 
to the high energy density of radiation in a SB region, however, a 
large fraction of the non-thermal energy is expected to be emitted 
as IC radiation at GeV - TeV energies. Since the cooling time for 
the IC process is so short in this environment, electrons with an en- 
ergy as low as ~ 20 (200) GeV are efficiently cooled for a PWNe 
age of 10 4 (10 3 ) years (Figure [T| right). 



3 EMISSION FROM THE PWN POPULATION 

The results obtained in the previous section suggest that the recent 
particle injection history is of most importance when studying the 
HE and VHE 7-ray emission from a population of PWNe in SB 
galaxies. In the following we will investigate how a population of 
PWNe might evolve in an environment similar to the one found 
in the SB nuclei of NGC 253 and M 82. The most important pulsar- 
related properties that impact the overall level of non-thermal emis- 
sion from PWNe in SB galaxies are (1) the number of pulsars that 
form PWNe, (2) the efficiency e with which these convert rotational 
spin-down energy to non-thermal particles, and (3) the distribution 
of pulsar birth periods, as these determine the time evolution of 
particle energy input. 



3.1 NGC 253 

Based on th e FIR luminosity in the SB nucleus of NGC 253, |Melo| 



et al. 



1 2002 1 derive a SFR of ~ 3.5 A/ Q yr. With a Salpeter ini- 



tial mass function and assuming that only stars with masses be- 
tween 8 A/0 and 20 A/0 undergo type-II SN explosions and hence 
form pulsars, we estimate a type-II SN rate in the SB nucle us of 
NGC 253 of 0.02 yr" 1 . This value is close to the estimate by 



gelbracht et al.| < |1998| > of 0.03 yr" 1 
the 0.08 yr" as estimated by 



En- 



but co nsiderably l ower t han 



Van Buren & Greenhouse 



(19941. In 



1 Note that |Atoyan & Ah aronian ( 1996] derive a factor of 2.24 for the case 
of a spherical emission region. 



the following we assume that all pulsars that are left behind by SN 
type-II explosions form a PWN and evolve as discussed in Sec- 
tion]^ A key parameter for estimating the flux from these systems 
is the typical birth period Po of pulsars, as the energy available for 
particle acceleration is oc P ~ 2 . Birth periods have been estimated 
for only a small number of objects, for example for the Crab pul- 
sar where Po w 16 ms. However, a relatively recent study based 
on the population of radio pulsars suggests a mean birth period of 
(Po) ~ 300 ms, with a spread of ap ~ 150 ms l |Faucher-7j iguere 
|& Kaspi|2006) . More recently, a model taking into account the pop- 
ulation of 7-ray pulsars as detected by the Fermi-LAT suggests a 
revision of the mean pulsar birth period down to (Po) ~ 50 ms 
(Watters & Romani 2011). The Crab may therefore represent a 
rather extreme case. In the following we assume a fixed value of 
Po = 35 ms and discuss the implications of this choice below. A 
conversion efficiency e = 0.5 is chosen as discussed above. The 
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Figure 1. Fraction of the total power injected in electrons that is radiated in different spectral bands for a PWN with a "Crab-like" birth period of Pq = 20 ms 
and a braking index n = 3, after 10 3 yrs (left) and 10 4 yrs (right). The PWN is modeled as evolving in a SB-like (solid line) and ISM-like environment (dashed 
line), with a constant magnetic field strength of (left) and 5(j,G (right). The ISM radiation fields considered are, the CMB, a cold dust component with 

T = 50 K and a starlight component with T = 5000 K. Both radiation fields are assumed to have = 1 eV cm - 3 . The SB radiation field has a temperature 
of T = 50 K and (7 p h = 2400 eVcm - 3 . Note the difference in current spin-down power E of a factor of ~ 50 between 10 3 yrs and 10 4 yrs. The dominance 
of IC emission in the SB environment at all times is clearly visible. 



magnetic field in a PWN is expected to decrease with time, how- 
ever, in the SB environment energy losses are likely dominated by 
IC emission at all times and the B-field evolution is not expected 
to result in a modified particle energy distribution. We therefore as- 
sume no magnetic field evolution in the nebulae and use a mean 
Bpwn field value of 25 /iG. As in the previous section we assume 
n = 3, and injection of particles with power-law index a = 2.0. 

The left-hand panel of Figure[2]shows a model spectral energy 
distribution (SED) for the total non-thermal emission of a popula- 
tion of PWNe in NGC253's SB at different evolutionary stages. 
This calculation is a two-zone, time-dependent model, where parti- 
cles are injected in each PWN according to Equation[T] The num- 
ber of PWN of a given age is calculated in 0.25 dex steps under 
the assumption of a constant pulsar birth rate as given above. In- 
jected particles cool via IC and synchrotron processes in the PWN 
environment. For objects older than 5 x 10 4 years, we assume that 
no more injection of particles in the PWN occurs, that all particles 
escape and continue to cool in the SB region of NGC253. Non- 
thermal emission generated by this population of escaped particles 
add to the diffuse emission of NGC253's SB region. As the SB 
nucleus is characterised by a much higher magnetic field and av- 
erage density (Bsb — 250 fiG, nn ~ 250 cm -3 ) synchrotron, 
Bremsstrahlung and Coulomb losses become important for parti- 
cles that have escaped the PWNe. The high thermal and cosmic- 
ray pressure in the SB region leads to the formation of a SB wind 
which effectively removes material from the centre of the galaxy. 
The height of the SB r egion of H ~ 60 pc and th e wind speed of 
^wind ~ 300kms _1 JZirakashvili & V61k|2006| implies an ad- 



vective loss time of r ac i = (H/2)/v vr i a d — 10 s yrs, after which 
particles are removed from the SB region. Hence, the SED for an 
age of 1.5 x 10 5 yrs may well be representative of the equilibrium 
spectrum produced by all PWNe over the past 150 kyr. The SED 
for an age of 5 x 10 4 yrs, on the other hand, shows the contribu- 
tion of non-thermal particles that have not yet escaped their PWNe. 



This SED therefore represents the contribution of individual PWN 
emission (rather than diffuse emission) to the overall non-thermal 
emission from the SB nucleus of NGC 253. 



The right-hand panel of Figure [2] illustrates how the different 
energy-loss mechanisms shape the resulting SEDs. In the employed 
model, IC emission dominates at photon energies above ~ 100 keV 
due to efficient cooling of non-thermal particles in the strong om- 
nipresent FIR radiation field. The peak at ~ 500 GeV in the y- 
ray spectrum is caused by particles that have left their low-density 
PWNe environments after 5 x 10 4 yrs and emit Bremsstrahlung in 
the dense gas of the SB region. The synchrotron emission spectrum 
is complex with a double-hump structure. Since the magnetic field 
in the SB region is so high, electrons with an energy of > 1 GeV 
that have left the PWNe predominantly lose energy via the syn- 
chrotron process and contribute to the radio emission. The second 
hump of the synchrotron spectrum in the X-ray domain is caused 
by the population of higher-energy particles that cool in the mag- 
netic field of young PWNe. This figure also compares the two-zone 
model described above to a simple one-zone model. In this ap- 
proach, particles are injected continuously rather than according 
to Equation [T] representing the semi-continuous nature of particle 
injection from the population of PWN as a whole. The environment 
in which particles cool in this calculation has been chosen to match 
the PWN environment for the target fields for processes that are 
important for the highest energy particles (B ~ 25 /j,G and SB ra- 
diation field), but using the mean SB density (ng ~ 250 cm -3 ) for 
the calculation of Bremsstrahlung, as particles are expected to es- 
cape their PWN on roughly the Bremsstrahlung cooling time in the 
SB. Although the IC and Bremsstrahlung parts of the 7-ray spectra 
are quite similar, the synchrotron components look very different. 
This is a direct consequence of the two-zone model, with the low 
(high) magnetic field strengths inside (outside) of PWNe. 
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Figure 2. Predicted SEDs for PWNe populations in the SB environment of NGC 253 at different evolutionary stages. Also shown is radio data from the Very 
Large Array (VLA, |Cari lli 1996 , the X-ray point is emission from the central source X34 as given by Pietsch et al. (2001) and the Fermi-LAT and H.E.S.S. 
data is from Abramows kTetafJ 2012) . Left: The integrated emission from individual PWNe is indicated by the long-dashed, grey curve. The equilibrium 
spectrum for all objects with ages less than 1.5 X 10 5 yrs is given by the solid black curve. The effect of particle confinement in the SB region over the age 
of the SB episode of 10 7 yrs is indicated by the short-dashed black curve. Right: Equilibrium spectrum of NGC 253 for the two-zone, time-dependent model 
(straight lines) and for the simple one-zone model (dashed lines), both of which are described in the main text. The different components of the emission are 
IC (red), Bremsstrahlung (green) and synchrotron (blue). 



3.2 M82 



In the following we apply the same two-zone model calculation to 
M 82 and compare it to multi-wavelength data from radio to VHE. 
The FIR emission from M 82's SB nucleus can be well described 
by a dominant dust component with a temperature o f ~ 50 K and 
a luminosity of 3.8 x 10 10 L (Colbert et al.|l999) . The inferred 
SFR is ~ 10 M Q yr _1 , considerably higher than in NGC 253, and 
implies a type-II SN rate of ~0.06yr _1 . Adopting the disc-like SB 
geometry of 150 pc radius and 60 pc height as used in Stric kland &| 
|Stevens| ( |2000| >, the higher intensity of the FIR radiation in M82's 
SB region also leads to a higher radiation field energy density of 
C^ph,M82 — 5000 eVcm~ 3 . The magnetic field adopted here is the 
same as in NGC253's SB: £>sb = 250/iG, consistent with esti- 
mates of 200/xG - 400/xG (Thompson eTaLl[2006) |Strickland &] 
Heckman 20071. For the average density we again assume a value 
of n H — 250 cm -3 . 



Figure|5]shows the result for the two-zone model as described 
above. Particles are injected according to Equation [T] and for the 
estimated SN rate of 0.06 yr _1 . In order not to overproduce the 
emission measured by VERITAS, the pulsar birth period has been 
reduced to 40 ms. All other pulsar-related parameters are kept as 
above. As before the time after which particles leave the PWNe into 
the surrounding SB environment is fixed to 5 x 10 4 years. Particles 
are again assumed to leave the SB region via energy-independent 
transport in the SB wind after 1.5 x 10 5 yrs. The calculated SEDs 
have very similar shape to the curves obtained for NGC 253, as 
expected due to the similar target radiation field, magnetic fields, 
and average particle densities in the SB regions. 




Energy (eV) 



Figure 3. As Figure [2] but for the emission from PWNe that evolve in 
the SB environment of M 82. Multi-wavelength data is from|Colbert et al" 
(1999| radio), |Strickland & Heckman] 42007| X-ray), |Ackermann et al. 
pOrll HE 7 rays) and |VERITAS Collaboration|p009l VHE 7 rays). 



4 DISCUSSION 

Figures [2] and [3] show that with a suitable choice of model param- 
eters the VHE 7-ray emission of NGC 253 and M 82 may be dom- 
inated by PWN as suggested by Mannhei m et al.| (2012) . In the 
following we discuss the plausibility of these parameter values. 

Magnetic and radiation fields. Figure [2] shows that young 
PWNe are expected to dominate the HE and VHE 7-ray emission. 
Magnetic field strengths in young objects such as the Crab Nebula 
are estimated to be ~ 100 — 400/nG (Atoyan & Aharonian|1996| 
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Zhang et al. 2008] |Volpi et al.|2008||Tanaka & Takahara|2010| >. The 
precise strengths of magnetic fields in these systems has little effect 
on the VHE fluxes unless the magnetic energy density U b exceeds 
[/ ph . Note PWNe that dominate in the Mannh eim et al.| ^2012) 
model have much lower magnetic fields ( de Jager et al. 2009). Our 
SB energy density estimates are based on a simplified geometry but 
are relatively robust, and exceed Ub as long as B < 500 /iG (or 
even higher for M 82). IC dominance over (at least most of) the 
early lifetime of pulsars in the SBs therefore seems likely. 

Particle Escape. In our model we assume that all particles es- 
cape the PWNe after 5 x 10 4 yrs. Smaller values of the escape time, 
and energy-dependent escape are certainly plausible (see e.g. |Hin-| 
|ton et al.|201 1| > but have modest impact on the resulting SED except 
to boost diffuse synchrotron emission at high radio frequencies and 
bremsstrahlung emission at ~ 1 GeV. The escape time of particles 
from the SB region must be less than ~ 1 Myr to avoid violating 
constraints on the diffuse radio emission. However, energy inde- 
pendent escape via the SB wind on much shorter timescales seems 
likely (Zirakashvili & Volk 20061. For even higher wind speeds 



(e.g. Strickland & Heckman 2009 1, particles would escape the SB 
region faster and the Bremsstrahlung and radio components would 
be lower. 

Birth Period. Average values of P (or strictly (< P 2 >) 1//2 ) 
of ~35 ms and ~40 ms are required to match the VHE fluxes of 
NGC253 and M82, respectively. Whilst these values are some- 
what lower than typical assumptions for pulsar populations, it is 
certainly possible that these values are representative of recently 
formed pulsars in these objects. Furthermore, a plausible increase 
of the injection efficiency e and/or a systematic underestimate in 
the SN rate, would increase these values into a more comfortable 
range. 

Injection spectrum. The assumed hard injection spectrum 
(a = 2) maximises the VHE flux. Softer spectra at injection are 
plausible, and would be consistent with the measured VHE spec- 
tra, but increase the required energy input dramatically. For exam- 
ple for a — 2.3 the required birth periods are ~15ms and fast 
particle escape (t aBC <C 10 5 years) is required for consistency with 
radio limits. As a consequence PWN cannot explain the emission 
detected with Fermi from these objects. 



5 CONCLUSIONS AND OUTLOOK 

Considering the very different evolution of particle populations in 
starburst environments to typical Galactic environments we never- 
theless reach a similar conclusion to that of Mannheim et al. (' 2012) : 
that PWN can plausibly contribute to the TeV emission of NGC 253 
and M 82. The Fermi emission cannot be explained in this frame- 
work on energetics grounds and the standard explanation of diffuse 
emission from strong interactions of accelerated cosmic rays re- 
mains the most likely possibility. As a cut-off in this component 
just before the TeV range seems unlikely, the observed VHE emis- 
sion may represent a superposition of a harder PWN component on 
the diffuse p-p emission. CTA (Actis et al.|201 lp will have the sen- 
sitivity to search for the spectral feature associated with PWN and 
place tight constraints on the pulsar population. 
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